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We present an exact diagonalisation study of bilayer quan- 
tum Hall systems at filling factor v = 2 in the spherical ge- 
ometry. We find the high-Zeeman-coupling phase boundary 
of the broken symmetry canted antiferromagnet is given ex- 
actly by previous Hartree-Fock mean-field theories, but that 
the state's stability at weak Zeeman coupling has been quali- 
tatively overestimated. In the absence of interlayer tunneling, 
degeneracies occur between total spin multiplets due to the 
Hamiltonian's invariance under independent spin-rotations in 
top and bottom two-dimensional electron layers. 

In the last decade there has been an increasing interest 
in quantum Hall ferromagnets 0. Most recently bilayer 
systems at a filling factor of two have become the object 
of intensive theoretical ||- || and experimental Q- 
fl2f research. The rich phenomenology of v = 2 bilayers 
mirrors a complex interplay between Coulomb interac- 
tions in the lowest Landau level, Fermi statistics, and 
the coupling of external fields to spin and layer degrees 
of freedom. 

Our current microscopic understanding of v = 2 bilayer 
quantum Hall ferromagnets is based on Hartree-Fock 
mean-field theory calculations and on a partially 

phenomenological effective spin Hamiltonian description 
H . Both approaches lead to the prediction of a novel bro- 
ken symmetry canted-antiferromagnet ground state with 
finite spin-suscpetibility which interpolates, as external 
field parameters are varied, between a fully spin-polarized 
state and a spin singlet state, both of which have charge 
and spin gaps and zero differential spin susceptibility. In 
this Letter we report on the first exact diagonalisation 
study of finite bilayer systems at filling factor v = 2. 
Our calculations support the predicted occurrence of a 
broken symmetry ground state, and indeed demonstrate 
that the Hartree-Fock result for its large Zeeman cou- 
pling phase boundary is exact. We find that the stability 
of the canted antiferromagnet state relative to the spin 
singlet state is overstated by the Hartree-Fock approx- 
imation, and estimate the correct position of the phase 
boundary. 

A bilayer system in a strong magnetic field is described 
in spherical geometry by the following Hamiltonian, 



H = H-ip + "Ho 



oul 



(1) 



where 7icoui represents the usual Coulomb interaction 
within and between layers, and the single-particle Hamil- 
tonian Hip is given by 



(2) 



A summation convention is understood for repeated 
greek indices, where /i, // G {+, ~ } run over the layer 
(or pseudospin) index, while a, a' £ {T^i} run over 
the z-projections of the electron spin; f, 3 are Pauli 
matrices for pseudospin and electron spin, respectively. 
m E {—N^/2, . . . , is the z-projection of the orbital 

angular momentum of each electron in the lowest Landau 
level, where iV^ is the number of flux quanta penetrating 
the sphere. The Hamiltonian contains bias voltage (A„), 
tunneling (At), and Zeeman coupling (A 2 ) terms. In the 
following we measure the interlayer separation d in units 
of the magnetic length / b = WKc/eB and all energies in 
units of the Coulomb energy scale e 2 /els- 
We first consider the case where all single-particle cou- 
pling constants vanish. The number of particles in each 
layer is then a good quantum number and, in the ground 
state, both layers have filling factor one. Moreover, be- 
cause the Coulomb interaction is spin-independent, the 
Hamiltonian is invariant under independent spin rota- 
tions in either layer. It follows that the total spin in 
either layer is a good quantum number. 
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FIG. 1. The low-lying spectrum of a double layer system 
at filling factor two with eight electrons. The energy levels are 
plotted as a function of the total orbital angular momentum 
L for different values of the total spin S. The upper diagram 
shows the case of vanishing single-particle couplings, while in 
the lower spectrum At is finite. For A t = 0, the degenerate 
spin multiplets in the ground state are exhibited explicitly . 



1 



For layer separation d — > oo it is known that the ground 
state of each isolated layer is a quantum Hall ferromagnet 
with 5 = N/4, where N is the total number of particles in 
the two layers. Since the Coulomb interaction within the 
layers is stronger than between them, we anticipate that 
this should remain true at any finite d. Indeed, this ex- 
pectation is confirmed numerically. The upper diagram 
of figure [l] shows the low-lying spectrum of a double layer 
system at filling factor two with eight electrons in the case 
of vanishing single-particle couplings. The ground state 
as described there consists of multiplets with total spin 
quantum number 5 varying from to 4, the total spin 
representations of a state with good spin quantum num- 
ber S = 2 in each layer. When an interlayer tunneling 
term is added to the Hamiltonian, the additional sym- 
metry responsible for these degeneracies is lost and only 
total spin is a good quantum number. 
The robustness of individual layer ground state spin 
quantum numbers against the effects of added interlayer 
interactions is not restricted to total filling factor two. In 
fact, in our finite-size numerical calculations, it holds at 
all combination of numbers of flux quanta and electrons 
we have checked. 
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FIG. 2. The same system as in figure [l], but with one elec- 
tron removed. The degenerate multiplets discussed in the text 
are shown in detail. 



For instance, consider the system discussed above, but 
with one electron either added or removed. (These two 
cases are equivalent by particle-hole symmetry.) In the 
ground state, one of the layers has a filling factor of one 
(a ground state with L = 0, 5 = 2), while the other one 
contains a hole and the low-lying states are organised in 
the well-known skyrmion branch having quantum num- 
bers L = 5 = 1/2,3/2 @. The spectrum of this double 
layer system is shown in the upper diagram of figure |^. 
The ground state carries quantum numbers L = 1/2 and 
5 = 3/2,5 = 5/2 and is the result of coupling the L = 0, 
5 = 2 ground state of the full layer to the L = S = 1/2 
ground state of the hole system. Since either layer can 
carry the hole, we get two copies of the degenerate multi- 
plets. The next higher degenerate group of multiplets is 



the result of coupling the ground state of the v = 1 layer 
to the L = 5 = 3/2 state of the hole layer, and again 
all multiplets are doubled. To confirm this interpreta- 
tion we have verified that that this degeneracy doubling 
is lifted by applying a bias voltage, as shown in the lower 
diagram of figure |2| These results demonstrate that in 
the absence of interlayer tunneling, bilayer states for v 
near 2 can be safely regarded as two single-layer quantum 
Hall ferromagnets whose coupling has only a quantitative 
significance, for example in changing the energies of the 
skyrmionic elementary charged excitations. 

With this established, we now focus on changes in 
the nature of the ground state at v = 2 as the single- 
particle coupling constants are varied. An arbitrarily 
small tunneling amplitude is enough to break the de- 
generacy among the different spin multiplets and make 
the spin singlet state with the most antiparallel electron 
spin structure and consequently the most parallel pseu- 
dospin structure the nondegenerate ground state. More 
precisely, the energy levels are found to be ordered by 
the total spin 5 with the difference between neighboring 
levels increasing with increasing 5, as shown in the lower 
diagram of figure [l]. We note that the lone maximally 
polarized 5 = N/2 multiplet is annihilated by the in- 
terlayer tunneling term in the Hamiltonian and has an 
eigenenergy which is independent of At . Turning on the 
Zeeman coupling A z at a given value of the tunneling 
A t does not change the eigenstates themselves, but only 
shifts their energies and breaks the degeneracy within 
each spin multiplet. 

These findings lead to the following scenario: With in- 
creasing A z a lower critical Zeeman coupling A^\A t ) 
is reached where the state with S z = S = 1 becomes 
the ground state, i. e. the system leaves the spin singlet 
phase. If A z is increased further, the ground state S z = 5 
quantum number increases monotonously until, at an up- 

(2) 

per critical value, A z (A t ) the fully spin-polarized state 
S z = 5 = N/2 is reached. 

Finite-size spectra are shown in figure |^ for several 
Zeeman couplings at At = 0.8. In the top diagram the 
system is in the spin singlet phase, while in the bottom 
the ground state is fully spin-polarized. In the narrow 
transition area A* 1 ' < A z < A z 2 ^ all low-lying states 
with S z — 0, . . . , N/2 have energies very close to each 
other. This holds also for a branch of states with angular 
momentum L > which appears to be separated by a 
gap from higher-lying parts of the spectrum. We identify 
this transition region with the canted antiferromagnetic 
phase first proposed by Zheng et al. on the basis of the 
unrestricted Hartree-Fock approximation a. 
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FIG. 3. Finite-size spectra for of a system of eight electrons 
as a function of the Zeeman coupling for tunneling amplitude 
At = 0.8. States with negative values of S z are omitted for 
clarity. 

The mean-field state, which breaks spin-rotational 
symmetry around the z axis can be constructed as a 
linear combination of the nearly degenerate exact eigen- 
states \S Z ) carrying quantum numbers L = and definite 
S z values. To analyse the spin structure perpendicular 
to the Zeeman axis we introduce spin operators for each 
layer separately, S^, /i G { + ,—}, where the total spin 
of the bilayer system is given by S = ^ M S^. Since the 



states \S Z ) belong to different multiplets, all matrix ele- 
ments of between them are zero, which means that 



(S Z \S+\S Z -1) 



-(S 2 



\s^\s z 



1) 



(3) 



Thus, the matrix elements of spin components perpendic- 
ular to the Zeeman axis have the same magnitude and 
opposite sign in opposite layers. Any wavepacket con- 
structed from these states will, like the mean-field state, 
have opposite transverse spin-polarization in the two lay- 
ers. 



Next let us consider the phase diagram, i. e. the 
functions Ai^Ai) and Ai^(At). In the thermody- 
namic limit these lines mark the boundaries between 
spin-polarized, spin-singlet, and canted antiferromagnet 
phases. They may be compared with the Hartree-Fock 
results obtained recently in the planar geometry for an 
infinite system To avoid unnecessary finite-size un- 
certainty, we have rederived the Hartree-Fock equations 
for the spherical geometry obtaining explicit expressions 
for finite systems. We find that the phase boundaries 
have the same form as for the infinite system jl: 



AW = y/A t (A t - 2F_) 
for A t > 2F_, otherwise a£ X) = 0, and 
Al 2 ) = 



A? 



F 2 — F_ 



(4) 



(5) 



In the present case, however, the exchange parameter F_ 
is size-dependent: 



F_ = 



N rl> 



1 



d B ^2% 



with 



1(a) = 



"1-1 5 



1(a) 



dx- 



l 



i d 2 



(6) 



(7) 



In figure Q we compare exact-diagonalization and 
Hartree-Fock finite-size phase boundaries for N = 
6, 8, 10, 12 electrons. Not unexpectedly, the region of 
the canted antiferromagnctic phase turns out to be much 
smaller than predicted by the Hartree-Fock theory. In- 
terestingly however, the two results for the phase bound- 

(2) 

ary A z (A t ) coincide within our numerical precision of 
at least 10 -12 . Since both the spin-polarized state and 
its elementary collective excitations, which go soft at the 
phase boundary, are described exactly jl4j by Hartree- 
Fock theory, this coincidence is not entirely unexpected. 
It does, however, very convincingly demonstrate that the 
spin-polarized to canted phase transition remains con- 
tinuous when quantum fluctuations are included. We 
conclude that the position of this phase boundary in 
the thermodynamic limit can be calculated exactly using 
the expressions given in reference B , or equivalently the 
— * oo limit of equations (^-(Q), adding finite-well- 
width and other sample-specific corrections as required. 
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A, A, 
FIG. 4. Phase diagrams at a layer separation of d = 1.0 
and zero bias voltage for different system sizes. The upper 
and lower critical Zeeman coupling as obtained from exact 
diagonalisation (ED) data are plotted as thick lines. The 
canted antiferromagnetic phase lies between them and sep- 
arates the spin-polarized phase (SP) from the spin singlet 
phase (SS). The finite-size Hartree-Fock (HF) results are given 
for comparision. The upper HF phase boundary coincides 
with the ED result (see text), while the lower phase boundary 
strongly overestimates the region of the canted antiferromag- 
netic ground state. 

The last panel (N = 12) shows in addition the the HF phase 
boundaries for the infinite system (where the upper one is as 
well exact) as dashed lines. The value of At below which the 
ED results for A^(At) decrease with system size is marked 
by a vertical line. 



The result for this upper boundary is in marked con- 
trast with our findings for the lower boundary. Here the 
Hartree-Fock approach leads to a canted antiferromag- 
netic region at all system sizes even in the absence of the 
Zeeman coupling. A similar conclusion was reached by 
Demler and Das Sarma using the effective spin theory 
|^| . In our finite-size exact diagonalization calculations, 
on the other hand, A* is always finite. For A t < 0.35, 
marked by a vertical line in the last panel of figure ||, 
Az decreases monotonously with system size, while it 
increases monotonously for larger values of A t . These 
observations guarantee that the canted antiferromag- 
netic phase predicted by Hartree-Fock theory is actually 
present in the infinite system (although clearly dimin- 
ished by quantum fluctuations). Moreover, our findings 
are consistent with a vanishing thermodynamic limit of 
Aa (At) and a nonzero spin susceptibility for A t < 0.35. 



Static correlation function calculations |15| at A 2 = are 
also consistent with a nonzero canted order parameter 
for A t < 0.35. These numerical results are thus consis- 
tent with a single intermediate phase which has both a 
finite spin-susceptibility and canted-antiferromagnet or- 
der. For At > 0.35, A^ increases slowly with system 
size, presumably saturating at a finite value smaller than 

(2) 

A z and leaving a narrow canted-antiferromagnet strip 



in the phase diagram. The maximum A t value at which 
the ordered state phase extends down to A z = is much 
smaller than the value At = 0.60, predicted by Hartree- 
Fock theory. These findings are illustrated along with the 
exact upper phase boundary in the right bottom panel 
of figure |J, which shows that the Hartree-Fock approxi- 
mation strongly overestimates the stability of the canted 
antiferromagnetic phase against the spin singlet phase in 
finite systems and as well in the thermodynamic limit. 
JS acknowledges support from the Deutsche Forschungs- 
gemeinschaft under grant SCHL 539/1-1. AHM was sup- 
ported from the National Science Foundation under grant 
DMR-9714055. 
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